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We present a unified account of damping of low-lying collective modes and of relaxation of tem- 
perature anisotropies in a trapped Bose gas in the collisionless regime. By means of variational 
techniques, we show that the relaxation times for the two situations are closely related to the 
simplest variational estimate of the viscous relaxation time. We derive rather precise theoretical 
expressions for the characteristic relaxation times, and compare our results with experiment. 
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Many experiments have been reported recently on re- 
laxation processes in clouds of bosons in magnetic traps. 
Among these are measurements of damping of collective 
modes in vapors of rubidium [QJ^] and sodium , and 
of relaxation of temperature anisotropies in cesium [^|,^| , 
rubidium and sodium vapors It is an impor- 

tant task for theory to understand these processes for a 
number of reasons. One is to interpret the experimen- 
tal data, and to use them to deduce atomic scattering 
properties. The second is theoretical, since these sys- 
tems are qualitatively different from most other atomic 
systems with which one is familiar, in that they are finite, 
and that the mean free path of an atom is, in the vast 
majority of experiments, comparable with or larger than 
the dimensions of the cloud. One is thus not in the hy- 
drodynamic regime, for which so much of the traditional 
kinetic theory of gases was developed. In attacking the 
problem of relaxation in trapped clouds of atoms, one 
has to confront a number of fundamental issues also of 
importance in other contexts, such as collective motion 
in nuclei Jl0|]. In this Letter we focus on the properties 
of clouds above the Bose-Einstein condensation temper- 
ature, T c . This work is of importance in its own right, as 
well as for understanding physical processes below T c . 

In the past, the theory of relaxation of temperature 
anisotropies in trapped clouds of alkali-atom vapors has 
been approached both analytically and numerically, 
using Monte-Carlo simulations (||3[] or other methods 
jn]], and we know of no calculations of the damping 
of modes in other than the hydrodynamic limit p^| , ^3| . 
Here we shall attack the problems analytically, the basic 
physical idea being to consider the motion of atoms in 
the absence of collisions, and then to treat the effect of 
collisions as a perturbation. 

We begin with some general considerations. The traps 
employed in experiments performed to date have an axis 
of symmetry, which we shall refer to as the z axis. Gen- 
erally the trapping potential, V(r), is harmonic, and we 
shall write it in the form 
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where uj z and wj_ are the frequencies associated with mo- 
tion of atoms in the trap along the axis of the trap, and 
perpendicular to it. The effect of the interaction en- 
ergy on the mode frequency is small, of relative order 
nUo/ '(ksT), where Uq — ATrh 2 a sc /m is the effective two- 
particle interaction. Here n is the particle density, a sc is 
the scattering length for atom-atom collisions, m is the 
mass of an atom, and ks is Boltzmann's constant. The 
quantity nUo/ (ksT) is typically of order one per cent, 
and we shall neglect it. In the absence of collisions, mo- 
tions along the symmetry axis and perpendicular to it 
are independent, with modes associated with motion in 
the z direction having frequencies which are multiples 
of lu z , and those associated with motion in the x and y 
directions having frequencies which are multiples of u)j_. 
Our approach is quite general, but since modes with fre- 
quency close to 2uj z and 2lu± have been investigated in 
detail experimentally, we shall focus on these. 

Quantum mechanically, the spectrum of single-particle 
levels is discrete. However, the distribution of levels will 
be spread out due to deviations of the trap potential 
from the purely harmonic form, and due to interactions 
between particles. We shall assume that for conditions 
of experimental relevance, with single-particle energies 
of order 10 2 oscillator quanta, the spreading is sufficient 
that the particles may be described semiclassically, in 
terms of a distribution function /(r, p), where p is the 
particle momentum. We now calculate the damping of 
collective modes above T c . 

Since the system is dilute, in the sense that the typi- 
cal particle separation is small compared with the atom- 
atom scattering length, the particle distribution function 
satisfies the Boltzmann equation, which is 
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where /[/] is the collision term. Let us now consider an 
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oscillation involving collective motion in the z direction. 
In the absence of collisions, the general solution of the 
Boltzmann equation having time dependence e~ 2iuJ ' t is 

Sf = zlgiKj), 



(3) 



where z± = p z ± imio z z are the classical equivalents of 
the raising and lowering operators in quantum mechan- 
ics, g is an arbitrary function, and Kj are constants of 
the motion. Since z± varies in time as e^ luJzt , it is ob- 
vious that Sf, Eq. (||) has the desired time dependence. 
One may ask whether more general forms are possible, 
based on, say However, = 2mE z where 

Ei = p 2 /(2m) + muifxf/2 is the energy associated with 
motion in the i direction. Thus this expression is still of 
the general form given by Eq. (||). 

In the absence of collisions there is enormous degener- 
acy among modes that have frequency 2lu z . This corre- 
sponds in the quantum mechanics of particles in a sim- 
ple harmonic oscillator potential to the many different 
ways of giving an excitation energy 2tuv z to the system. 
This degeneracy is broken by collisions between parti- 
cles, and we now wish to find the mode which has the 
longest lifetime, because we anticipate that this will be 
the one of experimental relevance. Collisions will also 
mix parts of the distribution with different time depen- 
dences in the absence of collisions, e.g. those for which 
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However, since we are interested in 



the collisionless regime, the collision rate is small com- 
pared with the oscillator frequency and this mixing will 
be small. It is convenient to introduce the notation 
$(r,p) = <S/(r,p)/[/°(l + /°)], where /° is the equilib- 
rium distribution function, in terms of which the Boltz- 
mann equation is given by 
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where the linearized collision integral is 
I{<P] = J dT i J d(T \ v - vi|($ + $i - $' - $i) x 

x/°/i (l + / ')(l + A '). (5) 

Here the momenta of the incoming particles in a collision 
are denoted by p and p', and those of the outgoing ones 
by pi and p' 1; and dr = dp/(2nh) 3 denotes the volume 
element in momentum space. We have introduced the 
differential cross section da, which in general depends on 
the relative velocity u = |v— vi| of the two incoming par- 
ticles as well as on the angle between the relative velocity 
before and after the collision. In this study we restrict 
ourselves to the case of s-wave scattering, since the ener- 
gies of the colliding particles are quite low; we shall also 
assume for the moment that the cross section is indepen- 
dent of the momentum, which is a good approximation 
for Rb and Na atoms, but later we shall examine the 
case of a momentum-dependent cross section, which is 
relevant for Cs. 



We now look for solutions of the Boltzmann equation 
having the form e —%«<>zt-rt _ These are j us t the eigenfunc- 
tions of the collision term, ^i, which satisfy the equation 



ri/°(i + /°)*i = /[*i], 



(0) 



where I\ is the corresponding eigenvalue. The most long- 
lived mode corresponds to the smallest value of I\ , which 
we denote by Tq. The task of finding Tq is most simply 
attacked by variational methods, since for any trial func- 
tion W it is clear that 



r < 
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which is our basic variational expression for the damp- 
ing rate. Here the brackets denote integration over both 
coordinate and momentum space. 

As a first application, we estimate the damping rate 
of the mode considered above. We employ the simplest 
possible trial function consistent with the fact that it 
must be of the form of Eq. (g) for modes with frequency 
2to z , and write W = z\. This trial function contains 
terms of the form z 1 , zp z , and p 2 z . As is clear from Eq. (g) , 
the collision integral gives zero when operating on the 
first two of these because of conservation in collisions of 
particle number and momentum, and, because of energy 
conservation, it also gives zero when operating on p 2 . 
Consequently we may write 

{{z* + ) 2 i[z 2 + ]) = (p 2 z i[ P 2 z }) = {(p 2 z - P 2 /mpi-p 2 m). (8) 

This is a particularly convenient form because precisely 
the same quantity occurs in the simplest variational cal- 
culation of the shear viscosity, r\, which is given by 

110, 

, vai - 1 [JMPi-P 2 m 2 f (i + f°)]\ (9) 

m 2 kBT J dr(p 2 — p 2 /3)I[p 2 — p 2 /3] 
Substitution of Eq. (||) into Eq. (0) leads to the result 

((p 2 z - P 2 /mp 2 -p 2 m) 



r n < 
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Let us begin by considering a non-degenerate gas. To 
bring out the physical content of our results we introduce 
the variational estimate for the viscous relaxation time. 
This is given by r 7))Var (r) = J?var/[n(r)fc B T], where n(r) 
is the local density. For the case of energy independent 
s-wave scattering one finds 04| 
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where a — 87rOg C is the total cross section, and v is 
the mean thermal velocity, v = (8/cbZ '/nm) 1 / 2 . From 
Eqs. (g) and @ we are led to the following expression 
for the damping 
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n(r) 
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dr n(r) 
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The factor of 1/6 reflects the fact that collisions damp di- 
rectly only the part of the distribution function varying 
as p 2 — p 2 /3, while all parts of the distribution contribute 
to the total energy associated with the oscillation. The 
average of the scattering rate that occurs here is propor- 
tional to / n 2 (r) dr/ J n(r) dr — n(0)/2 3 / 2 , where n(0) is 
the density at the center of the cloud, and therefore we 
find 



Vl-Vl = (87r/15)VV(r 2 2 +Kr 2 ) andp 2 +p 2 -2p 2 /3 = 
(2/3) 1 / 2 (8tt/15) 1 / V*2°) where Y™ are the spherical har- 
monics, we find that Ta(m z = 2) = 3ro(m z = 0), and 
thus, 



r (ni 2 = 2) < -^n ci (0)<7v 



(16) 



r n < 
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(13) 



One can systematically improve the estimates of mode 
damping by using more general trial functions. As in the 
case of transport properties of classical gases, one finds 
that the simplest variational calculation is extremely 
good, and yields results which are within a few per cent 
of the exact answer. Details of these calculations will be 
reported elsewhere. 

Including effects of degeneracy to leading order in 
T c /T, we find 



We now compare our results with experiment. The 
damping of oscillations along the axis of the trap has 
been measured by the MIT group ||, and at T = 2T C 
they found a damping time of about 80 ms. The fre- 
quency, v — lu/2tt, of the mode is « 35 Hz, somewhat 
below the frequency vc — ~ 38 Hz expected in 
the collisionless regime, but significantly above the value 
vh ~ y/12/5 v z « 29 Hz in the hydrodynamic regime 
pH . This is clear evidence that the system is in an in- 
termediate regime between the collisionless and hydrody- 
namic limits, but somewhat closer to the collisionless one. 
For N = 5 x 10 7 at T = 2T C and for cj ± w 2vr x 250 
Hz, lo z Ri 2tt x 19 Hz, classically the central density is 
- 2.9 x 10 13 cm" 3 . Taking a sc = 27.5 A |7), we esti- 
mate that r. 
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48 ms in the classical limit, or sa 46 
ms if the effects of degeneracy are included. We can 
allow in an approximate fashion for the fact that the 
system is not in the collisionless limit by using the sim- 
ple interpolation formula, Eqs. (21) and (22) of Ref. Jl2[ , 
which indicates that the damping rate is approximately 
(y — vn)j(yc — vh) times the value in the collisionless 
limit. This gives a theoretical estimate for the damping 
time of ~ 70 ms, in good agreement with the experimen- 
tal value of ~ 80 ms. 

In a more recent experiment of the MIT group [Q, 
axial oscillations were measured to have a damping time 
of w 50 ms at T w T c and a frequency v sa 1.75 z^ z . In 
this experiment N w 8 x 10 7 at T m T c , luj_ ps 2tt x 230 
Hz, and uj z w 2tt x 17 Hz. Thus n cl (0) w 6.3 x 10 13 
cm" 3 and T c w 1.9 fiK, so Tg^m = 0) w 30 ms for 
classical statistics. The effects of degeneracy are more 
pronounced in this case, and if we take only the leading 
term, as in Eq. (|l4|), we find a damping time of ~ 24 ms. 
This experiment is slightly closer to the hydrodynamic 
limit than to the collisionless one, since the factor (y — 
vn)l(yc~ v n) is ~ 0.44, and the damping time estimated 
as in the previous case is ~ 55 ms, again in excellent 
agreement with the experimental value, ~ 50 ms. 

We now turn to measurements of modes involving mo- 
tion perpendicular to the axis of the trap. The m z = 
and m z = 2 modes have both been investigated in ex- 
periments by the JILA group Q|. This experiment is 
in the collisionless regime, since the product of the fre- 
quency lu of the mode and the scattering time r sc is ~ 20, 
and the mean free path is at least four times the typ- 
ical dimension of the cloud. Taking N = 8 x 10 4 at 
T = 1.1 T c , wiw27tx 129 Hz, and uj z w 2tt%/8 x 129 Hz, 
we find that classically the central density is ~ 3.7 x 10 13 
Concerning the m z — 2 mode, the part of ^ that con- cm" 3 . Using a scattering length a sc = 53 A, we find from 



Here n c i(0) = Nuj\<jj z [m/2Trk B T] z / 2 , with N being the 
total number of atoms, is the central density calculated 
neglecting effects of degeneracy, and Ci a ) is t ne Riemann 
zeta function. The coefficient 3£(3)/16 is approximately 
equal to 0.225. 

We now consider oscillations in the x and y directions. 
An oscillation in the x direction alone is an cigcnstate in 
the absence of collisions, but when collisions are included, 
it will be strongly coupled to the degenerate oscillation 
in the y direction. It is clear from the rotational symme- 
try of the problem that the appropriate eigenstates are 
ones having simple properties under rotation about the z 
axis. The simplest trial functions are thus $ = x\ + y\ 
(corresponding to magnetic quantum number m z = 0) 
and 5* = x\ — y\ (m z = 2), where x± and y± are defined 
by analogy with z±. 

Starting with the m z = mode, we see from Eq. (^) 
that this choice for $ gives the same answer for the colli- 
sion integral as before, since the only non- vanishing part 
comes from p 2 , +p 2 , which is equal to p 2 —p 2 z . The denom- 
inator of Eq. (0), though, is twice as large in this case. 
Since the correction due to the quantum degeneracy does 
not depend on the symmetry of the mode, the damping 
rate is 



T (m z = 0) < ^n ci (0)av 
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tributes to the collision integral is p 2 



p 2 , and since 



Eq. (|T^) that for the m z = mode the damping time is 
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r ( ^" 1 (m 2 = 0) « 100 ms, while for the m z = 2 mode it 
is T' 1 (m z — 2) w 33 ms. Our results agree in order of 
magnitude with the experimental results, but a more de- 
tailed comparison between theory and experiment should 
await a reappraisal of the experimental data. 

In experiments with cooled atoms, one often encoun- 
ters situations where the degrees of freedom associated 
with the motion of atoms in the various directions in the 
trap are out of equilibrium with each other . Crudely 
speaking, one may talk of anisotropies of the temperature 
of particles. We now show how temperature relaxation 
may be treated in our approach, and conclude that again 
the characteristic relaxation time is simply related to the 
viscous relaxation time that determines the damping of 
oscillations. If one changes the temperature associated 
with particle motion in the z direction by an amount ST 
and that associated with motion in the transverse direc- 
tion by an amount —ST/ 2, thereby conserving the aver- 
age particle energy, the change in the particle distribution 
function is proportional to f°(l + f°)[E z ~(E x +E y )/2} cx 
— (x + x^ + y + y_)/2. The damping of such a dis- 
turbance may be treated in just the same way as for the 
case of oscillations, and one finds that the temperature 
relaxation rate, IV, is simply that given by Eq. (|Tl]), for a 
density n = n c \ (0)/2 3 / 2 , multiplied by 1/2. The physical 
reason for this factor of 1/2, is that both kinetic and po- 
tential energies have to be dissipated through collisions, 
whereas collisions degrade directly only the anisotropies 
associated with the kinetic energy. Therefore, we find for 
IV up to leading order in T c /T 



I> < ^n ci (0)(Tv 
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The classical version of the above equation can be ex- 
tracted from Ref. ||. The number 5 which appears in 
the coefficient of the above equation should be compared 
with the number 5.4 that results from a Monte Carlo 
analysis ||. 

Our results may readily be generalized to more general 
forms of the scattering cross section, a situation relevant 
for Cs which has a two-body resonance close to zero en- 
ergy ||. Assuming a total s-wave cross section cr(u), one 
finds in the classical limit a temperature relaxation rate 
given by 



r T < 



l 
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n c \(0)a(u)u, 



where 



a(u)u 



J °° duu e exp(—mu 2 /4:kBT)a(u)u 
Jp 00 duu 6 exp(— mu 2 / ^ksT) 



(18) 



(19) 



Under the experimental conditions Q it is a rather good 
approximation to assume that the scattering is resonant, 
with a — 87r/fc 2 where k = mu/(2h) is the relative mo- 
mentum. We then obtain the result 



r T < 



128 n c i(0)^ 2 



15 



(20) 



The coefficient 15 should be compared with the value 
2 x 10.7 w 21 obtained from a Monte Carlo simulation 
1- 
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